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Abstract 

We compute the electromagnetic mass differences of mesons containing a single heavy quark in terms 
of measurable data using QCD-based arguments in heavy-quark effective theory. We derive an un¬ 
subtracted dispersion relation that shows that the mass differences are calculable in terms of the 
properties of the lowest-lying physical intermediate states. We then consider the problem in the large- 
N limit, where N is the number of QCD colors. In this limit, we can write a kind of double-dispersion 
relation for the amplitude required to determine the electromagnetic mass difference. We use this to 
derive analogs of the Weinberg sum rules for heavy meson matrix elements valid to leading order in 
\/N and to 0(1 /toq) in the heavy quark expansion. In order to obtain our final result, we assume 
that the electromagnetic mass differences and sum rules are dominated by the lowest-lying states in 
analogy with the situation for the 7r+-7r° mass difference. Despite the fact that some of the matrix 
elements appearing in our final result have not yet been accurately measured, we can obtain useful 
estimates: for example, we obtain {M^+ — ~ -1-1.8 MeV. We argue that our results are 

accurate to about 30%. 


* e-mail: luty@ctp.mit.edu 
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1. Introduction 


The computation of the electromagnetic mass differences of hadrons is one of the classic problems 
in strong interaction physics (see e.g. [1][2][3]; for reviews, see [4] [5] [6]). At one time, it was believed 
that electromagnetism was the only source of isospin breaking, so that computations of electromagnetic 
mass differences could be compared directly with experiment. With the advent of QCD, it is now 
understood that isospin breaking arises both from electromagnetism and from differences in the u and 
d quark current masses. The modern motivation for computing electromagnetic mass differences is to 
disentangle the electromagnetic contributions to isospin-violating mass differences from those of the 
current quark masses in order to obtain information about the current quark masses.* 

In this paper, we will compute the O(e^) electromagnetic mass differences of lowest-lying mesons 
with quantum numbers Qu and Qd (denoted here by P-^ = u or d). We work to 0(1 /toq) in the 

heavy-quark expansion and to leading order in the 1/7V expansion, where N is the number of QCD 
colors. We give a detailed exposition of the formalism used and give useful estimates of the electro¬ 
magnetic mass differences.^ A detailed comparison to experiment and the extraction of information 
about the light quark masses is carried out in a separate paper [8]. 

We follow as closely as possible the method of the classic calculation of the 7r+-7r° mass difference 
[2] [6]. The basic strategy is to write the O(e^) self-energy of the meson in terms of the forward Compton 
amplitude T, and then write a dispersion relation for T expressing it in terms of measurable data. 
In the case of the 7r+-7r° mass difference, one can use the fact that the pions are pseudo-Nambu- 
Goldstone bosons to write T in terms of vacuum current correlation functions depending only on a 
single kinematic invariant. Dispersion relations then relate these to data measured in e“'"e“ ^ hadrons 
and r decays. For heavy mesons we must work directly with T, which depends on two kinematic 
invariants. This makes the computation of the heavy meson electromagnetic mass differences more 
complicated than the 7r+-7r° mass difference. 

We begin in section 2 by writing an unsubtracted dispersion relation that shows that T is deter¬ 
mined by the properties of low-lying meson states and the lowest-lying excitations of the heavy-light 
quark system. Unfortunately, this dispersion relation cannot be used directly to compute T; one rea¬ 
son is that it depends on structure functions for timelike photon momenta, which are not measurable 
in practice. In section 3, we consider the problem in the combined heavy-quark and large-limits. 
In the large-limit, we can write T in terms of heavy-meson form factors and meson scattering am¬ 
plitudes. We derive sum rules relating the properties of the masses and matrix elements of the states 
that determine T. (These are exact analogs of the Weinberg sum rules for vector- and axial-vector 
correlation functions in the vacuum [9].) By assuming that these sum rules are approximately satu¬ 
rated by the lowest-lying states, we can compute the heavy meson electromagnetic mass differences 
in terms of heavy meson form factors. This is analogous to the successful classic calculation of the 
7r+-7r° mass difference. In section 4, we consider the l/mg corrections to these results, and section 
5 summarizes our results and gives our conclusions. This is a rather long paper, but the reader can 
get the main ideas by reading sections 2 and section 3 through subsection 3.1, followed by section 5 

* The 7r+-7r° mass difference is a special case, since quark masses do not contribute to it at leading 
order. Therefore, the electromagnetic contribution is expected to dominate, and it can be compared 
directly with experiment. 

f The electromagnetic mass differences of heavy mesons have also been considered in ref. [7], but 
the methods used are not based on a systematic approximation of QCD. 
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(which contains a summary of the main results). The sections omitted in this way consist mainly of 
repeated application of the ideas in the first part of the paper. 

Our final result is similar to the prescription used long ago to compute baryon electromagnetic 
mass differences [1], but we emphasize that this prescription was never put on a firm foundation [4] [5]. 
In fact, it is ironic that we are not able to extend our results to baryons because the large-limit is 
more complicated for baryons. 


2. Technical Preliminaries 

In this section, we derive some technical results that provide the foundation for the rest of the 
paper. Although we focus on the electromagnetic mass differences of heavy mesons in this paper, 
most of the formalism in this section can be applied to any type of hadron. 


2.1. Renormalization and Finiteness 


The O(e^) electromagnetic contribution to the mass difference of the pseudoscalar mesons P-^ 
{(. = u, d) with flavor quantum numbers Q£ is given by 


AM = Mu-M-^ = — 
“ d 2 


r d'^q AT{p,q) 

J (27r)4 g2 _|_ jQ_|_ 

- {P^ip)\dC\P^{p)) + {P-^ip)\SC\P-^{p)), 


( 2 . 1 ) 


where SC is the O(e^) counterterm in the underlying lagrangian required to render the result finite 
(see below); M is the mass of the heavy mesons in the absence of electromagnetic interactions; we 
have defined 

AT = Tu- T-^, (2.2) 

where 

T-,{p,q) = i Jd‘^xP’^ ^P-,ip)\TJ^^{0)J^{x)\P-,ip)). (2.3) 

Here, is the electromagnetic current, so T is a trace over the forward Compton amplitude for scat¬ 
tering of photons from heavy mesons. For states containing heavy mesons, we use the normalization 


{P{p)\P{p')) = {2TTf5^{p-p') 


(2.4) 


appropriate for the heavy particle effective theory. Similar equations hold for the lowest-lying heavy- 
light vector mesons P* that are related to P by heavy-quark symmetry. 

We now consider the counterterm SC in eq. (2.1). Because we will use the heavy-quark expansion, 
the underlying lagrangian is that of heavy quark effective theory [10]. To ©(l/mg) in the heavy-quark 
expansion, the lagrangian is 

C = - j tr(G^"G^,) + qilpq + Qiv ■ DQ 

1 _ _ h - (2-^) 

- - QD'^Q +^F^''Qa^,Q + ■■■. 

2mQ Amq Amq 

Here, is the electromagnetic field strength, is the gluon field strength and v is the 4-velocity 
of the heavy meson. We set the light current quark masses to zero for our computation, since we are 
not interested in 0{e^mu,d,s) effects. The coefficient of the term QD^Q is fixed by reparameterization 
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invariance [11]; the coefficients a and b can be fixed by matching in QCD perturbation theory and are 
unity to leading order in as{mQ) [12]. 

All divergences that appear when computing physical quantities to 0(l/mQ) can be absorbed 
into counterterms of the form appearing in eq. (2.5). This severely restricts the ultraviolet behavior of 
the matrix elements that appear in eq. (2.1), since only isospin-violating counterterms can contribute 
to the mass difference. Since we are neglecting light current quark masses, the only isospin-violating 
term in the lagrangian is the light-quark kinetic term qilpq. (It violates isospin because the light 
quarks have different charges.) However, matrix elements of qilpq between physical states vanish by 
the equations of motion. Thus, for the Qu-Qd mass difference, there is no counterterm contribution, 
and the integral on the right-hand side of eq. (2.1) must converge. (A similar argument in the context 
of full QCD is given in ref. [13].) In fact, since the quark charges are arbitrary parameters, the 
contributions from ATqq and ^Tgq must individually converge. 

The electromagnetic mass differences get contributions from isospin-violating counterterms at 
order Ijvr^Q, such as 

dL = ^F^’'q-/^D,q. ( 2 . 6 ) 

The integral for AM is convergent in full QCD [13], so this simply means that the 0(1 /toq) contri¬ 

bution to the electromagnetic mass difference is sensitive to momenta up to ~ mg, but the integrand 
falls off for momenta above mg fast enough so that the integral converges. 

It is convenient to rewrite AT in a useful form using some elementary isospin group theory. We 
first split the the electromagnetic current into heavy-quark and light-quark components 

J^ = J^ + J^. (2.7) 

The heavy-quark current is 

Jq = QqWQ, (2.8) 

where Qc = |, Qt = —| ■ The light-quark current can be decomposed into isospin 0 and 1 components 


■Jq — JqO + Jql > 


(2.9) 

where 

Jfj, Q'U Q'd _ n tU Q'U Q'd — n 1 

■JqO = 2 ^ 2 ^ - \ 


(2.10) 

and Q„ = 1, = —1. We can then write 



AT = ATqq + ATqg 


(2.11) 

in terms of the state (T = (Tgs alone: 



ATqq{p,q) = 2z Jd^xe^<^^P{p)\TJ^,iO)Jq,^{x)\P{p)) + [q- 


(2.12) 

ATQq{p,q) = 2z jd^xP‘>--{P{p)\TJ^{0)Jqi^{x)\P{p)) + {q- 

^-q). 

(2.13) 


We denote the corresponding contributions to AM by AMqq and AMqq, respectively. 
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2.2. Dispersion Relations and Low-energy Dominance 


Eq. (2.1) reduces the problem of computing the electromagnetic mass differences to the problem 
of determining the forward Compton amplitude AT defined in eq. (2.2). We address this problem by 
writing a dispersion relation that expresses AT in terms of measurable data. AT is a function of two 
kinematic invariants, which we take to be and v = va the frame where p = (Mp,0). We can 
then write a spectral representation for AT by inserting a complete set of states into eq. (2.3). The 
result is 


AT(p, 9^) 


2 ^^ , n'lmATin',q^) 

Jo + *0+ 


(2.14) 


which can be thought of as a fixed-tf^ dispersion relation. Dispersion relations such as this in general 
require subtraction because of the ultraviolet behavior of products of currents. However, due to 
isospin invariance, there is no counterterm that modifies the product of currents appearing in AT, so 
the operation of inserting a complete set of states that was used to derive eq. (2.14) is valid without 
modification.* 


ImAT is determined by physical matrix elements. In the mg ^oo limit, the relation is particu¬ 
larly simple: 


ImAT(p,g^) = 2t:^5(v - A„)(P(p)|J^(0)|n(g))(n((f)|Ji^(0)|P(p)) -bh.c., (2.15) 

n 

where M A„ is the mass of the state n, and Jq = = Jqi- The sum over n is 

over a complete set of states with fixed 3-momentum q. This in principle solves the problem of 
determining AT in terms of physical data, since ImAT(p, g^) is determined by structure functions 
for inelastic scattering off the heavy meson. However, eqs. (2.14) and (2.15) by themselves do not 
give a practical method of determining AT; for example, experimentally inaccessible data for timelike 
photon momenta (“timelike structure functions”) are required on the right-hand side of eq. (2.14). 
We will see in the next section that progress can still be made in the context of the large-A limit of 
QCD. 

The main use of the dispersion relation eq. (2.14) for our purposes is that it makes manifest 
that the electromagnetic mass differences are insensitive to the properties of very heavy states. We 
have argued that the integral that determines AM converges, so it is dominated by A.T{v,q^) with 
q"^ -^QCD- Now consider the contributions of an intermediate state n to AT in this kinematic 
region: the delta function in eq. (2.15) means that a state with mass M A„ contributes at ly' = A„ 
in the right-hand side of eq. (2.14). Therefore, the fact that the integral over ly' in eq. (2.14) converges 
means that the contribution of states with A„ ^ Aqcd is suppressed. 

The dispersion relation eq. (2.14) also shows that we can continue the integral in eq. (2.1) to 
Euclidean momenta. (Specifically, we write q = (iteo;?); so that d'^q = id'^qE and q^ = —q%.) We 

* In more detail: the operators are defined by functional differentiation with respect to ap¬ 
propriate sources S'g ^ added to the lagrangian. All divergences in this lagrangian can be canceled by 
counterterms depending on the sources with dimension 4 or less. Because there is no counterterm of 
the form 5'o5'i allowed by isospin invariance, there is no short-distance singularity in the product of 
currents appearing in AT, and the manipulations used to derive eq. (2.14) are valid. It is interesting 
to note that it was guessed that eq. (2.14) did not require subtraction before the advent of QCD; see 
for example ref. [4]. 
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also note that because the integrand is integrated over d'^qs, we can angularly average it in Euclidean 
momentum space without changing the result. We then have 

{AT){q%) 


AM = — 
2 


(27r)4 


4 


(2.16) 


where (AT) is the angular average of AT. Eq. (2.16) shows that we need not determine the full 
kinematic behavior of AT. For example, (AT) can be shown to be insensitive to the scaling behavior 
in the deep-inelastic region, in agreement with the fact that infinitely many states contribute to scaling 
behavior [5]. 


This problem of computing hadron electromagnetic mass differences has a long and involved 
history [4] [5], and we will not give a detailed discussion of the literature. We do wish to remark that 
one common approach to the problem of baryon electromagnetic mass differences has been to write 
a fixed-g^ (as opposed to fixed-(f^) dispersion relation and attempt to express the electromagnetic 
mass differences in terms of measured spacelike structure functions [3] [4] [5]. However, there are good 
reasons to believe that the fixed-g^ dispersion relation requires subtraction. (The necessity of the 
subtraction is related to the small x behavior of the structure functions; see e.g. ref. [5].) Since the 
electromagnetic splittings cannot be determined without knowing the subtraction constant, we do not 
follow this approach in this paper. We will, however, make some brief comments about the relation 
between this approach and ours in the next section. 


3. AM in the Large-A and Heavy-quark Limit 

In this section, we consider the computation of AM in the combined large-A and heavy-quark 
limit. We will consider 0(l/mQ) corrections in the next section, but we will not attempt to go beyond 
leading order in the 1/A expansion in this paper. 

For large A, QCD reduces to a weakly-coupled field theory with infinitely many meson fields 
whose interactions are polynomials in momenta [14]. The transition amplitudes can be expanded 
systematically in powers of 1/A, and the leading term in this expansion corresponds to keeping only 
tree graphs in the mesonic theory. The meson graphs that contribute to AT are shown in fig. 1. 

The representation of AT as a sum of tree graphs gives a kind of double-dispersion relation that 
will allow us to determine the heavy-meson electromagnetic mass differences in terms of measurable 
matrix elements. For example, it is clear that the graphs with intermediate heavy-meson lines (the first 
graph in fig. 1) are related to heavy meson form factors. However, we must also obtain information 
about the remaining graphs, which are not obviously related to form factors. We will do this by 
imposing consistency between the hadronic theory and properties of QCD and heavy-quark effective 
theory, such as those discussed in the previous section. These consistency conditions are expressed in 
terms of sum rules relating the (infinitely many) couplings of the states that appear in the graphs. If 
we assume that the sum rules are saturated by the lowest-lying states (consistent with the low-energy 
dominance proved above), we will find that the unknown contributions from the graphs in fig. 1 with 
no heavy-meson intermediate states give a numerically negligible contribution to the electromagnetic 
mass difference. 

For baryons, meson loops are not suppressed in the large-A limit (see [14][15]), so the arguments 
in this section cannot be simply extended to this case. Recent progress in the 1/A expansion for 
baryons [15] may be relevant to overcoming this difficulty. 

We now consider the contributions AMqq and AMg^ (defined by eqs. (2.12) and (2.13)) in turn. 
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3.1. AM, 


qg 


Consider AMgg from the point of view of the large-mesonic field theory. The sum of graphs 
that gives the amplitude ATgg can be written 


Arri f 2 2\ o AnWnil^"^, q^) 2 2\ 

AT„(. ,, ) = +2C(,. ,2 ), 


(3.1) 


where the first term is the contribution from the first graph in fig. 1, and the second term is the sum of 
the remaining graphs. Because ATgg is even under q —q, we consider it a function of rather than 
12 . The sum on n in the first term runs over the (infinitely many) heavy mesons with mass M -|- A„ 
that appear in the intermediate-state propagators in the first graph of fig. 1. The pole structure in 
this term arises by combining the heavy-meson propagators l/(±i^ — A„ -|- t0-|-). The vertices in the 
large-A limit are polynomials, and so the only non-polynomial dependence of Wn and C in eq. (3.1) 
comes from the meson propagators in figs. 1 and 2, which give poles in q^. Therefore, Wn and C are 
polynomial in 

We can give a physical interpretation of Wn{i2‘^, q^) by noting that for q^ spacelike, the only real 
intermediate states that can appear in ATgg in the large-A limit are single heavy meson states. From 
eq. (3.1) we compute 

lwLATgg{v'^,q^) = TT^Wn{q^)5{v - A„), (3.2) 

n 

where we have defined Wn{q^) = Wn{A^^q^). Substituting this into the spectral representation 
eq. (2.15), we obtain 


Wn{q^) = 2(P(p)| J^o(0)|n(g))(n(g)| Jgi^(0)|P(p)) -b h.c. (3.3) 


Therefore, Wn{q^) is proportional to on-shell heavy-meson form factors. (We can also see this directly 
by looking at diagrams.) 

Because the functions Wn are polynomial in we can expand Wn around = A^ to obtain 


^Tgg{v\q^) 


2 ^ AnWnjq^) 

^ _|_ ^Q_|_ 

rt. ■- 


Dn{^\q^) 


+ 2C{y\q% 


(3.4) 


where Dn{v‘^,q^) is polynomial in . The terms in eq. (3.4) proportional to Wn{q^) depend only on 
form factors, and we will refer to the contribution of these terms as the “form-factor” contribution. 
We refer to the remaining terms as “contact” contributions. (The first term in the sum over n is 
exactly what one would write for the right-hand side of an unsubtracted fixed-g^ dispersion relation. 
As mentioned in section 2, the fixed-g^ dispersion relation requires subtraction, and so we expect the 
two terms in the sum on n to separately diverge, while the total sum is finite. For this reason, we do 
not combine Dn with C.) 

The electromagnetic mass difference can therefore be written 


AMgg = 


fd^E 1 

J (27r)4 g| 


f AnWn{ — qE) 

\n + A2 


{Dn){ql] 


{C){ql) 


(3.5) 


where we have continued the integral to Euclidean momentum space. (Recall that (•) denotes the 
4-dimensional angular average.) 
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In the large-A^ limit, the functions C{v^^q^) and Dn{v^,q^) have the form 


C{v\q^), D^iiy^q^) 




polynomial in q^ 
polynomial in q^ ’ 


(3.6) 


where the polynomials in the denominator come from the vector-meson propagators in fig. 1. There¬ 
fore, after angular averaging in Euclidean space, we have 


(C')('7l) — ^2 


a 


( 0 ) 


/^(l) „2 

^rs HE 


iq% + ml^){q% + mi^) 




(3.7) 


{Dn){qE) — ^2 


D 


( 0 ) 


J-^nrsH E 


(<?: 


‘Or))?!; 




Pi?s\<ll) 


(3.8) 


where the sum over r, s runs over the vector mesons that appear in the propagators in fig. 1; these 
have isospin 0 and 1, respectively. (The contributions from terms that do not have two vector-meson 
poles can be written in this form by choosing the polynomials in the numerator to cancel the pole 
factors appearing in the denominator.) In eqs. (3.7) and (3.8), ..., are constants, and p(^'> 

and are polynomials in q^. 


The unknown quantities on the right-hand sides of eqs. (3.7) and (3.8) determine the contact 
contribution to the electromagnetic mass difference. Our strategy will be to restrict these quantities by 
imposing various consistency conditions arising from the definition of ^Tqq in QCD. These consistency 
conditions will involve the contact contribution alone because the form factor contribution will satisfy 
the consistency conditions by itself. 


We first note that the integral in eq. (3.5) must converge (as shown in section 2). It is expected 
on very general grounds that the form factors contributing to Wn{q^) fall off at large q^ (see below), 
so that for each n the form-factor contribution in eq. (3.5) is convergent. On the other hand, for 
general values of the coefficients in eqs. (3.7) and (3.8), the contact contribution will diverge in the 
ultraviolet, so the condition that the integral for AM is finite will constrain the contact contribution. 
If we impose a cutoff A on the photon momentum (note that this is gauge invariant), the divergences 
then have the simple forms {k = 1,2,...) and log A^. Demanding that the coefficients of these 
divergences separately vanish gives the “ultraviolet” sum rules 


E Pr? (?!) + E E ('?!) = 0’ (3.9) 

r,s n r,s 


E^-+EE^"- = o- (3.10) 

r.s n r,s 

Note that eq. (3.9) simply means that the polynomials P^^'> and pP'> do not contribute to AM. To 
derive these sum rules, it is crucial that the meson vertices are polynomial in momenta, so that there 
is a maximum possible power divergence A^^. If arbitrarily high powers of A^ were possible, then we 
would not be able to unambiguously separate the logarithmic and power divergences, since an infinite 
series in A^ may behave asymptotically like InA^.* 

* If the vertices that contribute to Dn are polynomials of higher degree as n becomes large, the 
integral in eq. (3.5) can give rise to arbitrarily high powers of A^. However, if we make some reasonable 
assumptions about the asymptotic properties of correlation functions in QCD, this can be shown not 
to occur. 
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Next, we note that as q—^0, the form factors simply measure the appropriate charge, so from 
eq. (3.3) we have 

Wniq^) ^{Ql - Ql)SnP a.sq—^0. (3.11) 

Substituting into eq. (3.4), we obtain 


ATggip,q)^27niQl-Ql)6iqo) + 2Y,DniO,0) + 2CiO,0) asg^O. (3.12) 

n 

(In the first term, we have used the identity 1/(^0 + *0+) + {q^—q) = —27rii5(<7o) to rewrite the 
P propagator.) The first term by itself is the correct result as q^O, as can be seen by writing a 
low-energy effective theory in which only the P and the P* appear. Taking the ultraviolet sum rule 
eq. (3.9) into account, we obtain the “infrared” sum rule 


E 


^( 0 ) 

L/rs 


EE: 


D. 


( 0 ) 




= 0 . 


(3.13) 


For the expert, we note that the sum rules in eqs. (3.10) and (3.13) are exact analogs of the 
Weinberg sum rules used in the calculation of the 7r+-7r° mass difference, in the sense that the 
Weinberg sum rules can be derived from identical considerations applied to the vacuum correlation 
function that appears in that calculation [6]. 

So far, all of the approximations we have made have been controlled, i.e. they become arbi¬ 
trarily accurate as some parameters of the underlying theory approach limiting values (N ^ 1 and 
TOQ Aqcd)- The result in eq. (3.5) is an infinite sum over one-particle states with infinitely many 
unknown parameters, even after the sum rules have been imposed. However, we know from the 
fixed-(f^ dispersion relation discussed in section 2 that intermediate states with large mass do not 
contribute significantly to AM. We do not know how many terms in eq. (3.5) are needed to get a 
good approximation to the sum, since there is no known small parameter controlling the convergence 
of the series. 

In order to make progress, we simply assume that the sum is well approximated by the first non¬ 
trivial term, i.e. that it is a good approximation to retaining only the mimimal set of intermediate 
states that gives a consistent description of the matrix elements that appear in the calculation. In 
particular, it is crucial that the matrix elements have the correct ultraviolet and infrared behavior to 
satisfy the sum rules derived above. 

Clearly, we must include both the P and the P* as intermediate states, since they become 
degenerate in the heavy-quark limit. To see what other states we must include, we consider the 
matrix elements 


{P{p)\J^o,imPiP + q)) = Fo,i{-q^)v^, (3.14) 

(^(p)k9\i(0)l^*(P + '?,e)) (3.15) 

which determine the form factor contribution if we include no heavy mesons other than the P and 
P*. In the appendix, we use the constituent counting rules [16] to show that the form factors F 
and G fall off as ~ l/g^ for large spacelike momentum transfer q. In the large-fV limit, the falloff 
of these form factors is due to the presence of single vector-meson intermediate states (see fig. 2). 
Since a single vector-meson pole can fall off at most as ~ l/(?^, we see that the asymptotic behavior 



requires cancelation between at least two different vector-meson intermediate states. The minimal set 
of vector-meson states we must consider therefore consists of the lightest isospin -1 vector mesons p 
and p' and the corresponding isospin-0 mesons u and lo' . From the behavior of the form factors at 
small and large q, we can write 


Qu ± Qd /3/2 

2 (l + 4/m2)(l + 4/m2,)’ 


(3.16) 

(3.17) 


where we have used the fact that mp = m^: and Wp' = rrit^i in the large-iV limit. (Numerically, 
rriuj — mp = 14 MeV.) Here, /3 is a strong-interaction matrix element that can be determined from 
the P* Pj decay rate. It is important to note that eqs. (3.16) and (3.17) are to be understood as 
valid only for the small q% that dominate in the integral that determines AM. 


We therefore truncate the intermediate states by keeping only the contributions from the P, P *, 
p, oj, p', and uj'. The contribution to AMqq from the form factor contribution to eq. (3.5) can be 
directly evaluated using eqs. (3.3) and (3.14)-(3.17). We have used the language of dispersion theory 
to make the connection to physical matrix elements explicit, but our Hnal result can be stated very 
simply in terms of Feynman graphs: we evaluate the one-loop electromagnetic self-energy for the 
heavy meson with momentum-dependent photon couplings as given in eqs. (3.16) and (3.17). The 
result is 

where 

a; = -^=0.53. (3.19) 

mp' 

(We use the symbol to denote statements that are valid only with the truncation of states 
discussed above.) Numerically, 

A^formfactor^-^0.38- 0.039 f- - —MeV. (3.20) 

qq GeV“^ J ’ 


1 -I- 3a: -I- 




(1 


a;)3 


8 (l-ba;)3 


(3.18) 


The value of (3 has not yet been measured directly. (There is a weak upper bound that results from 
using SU{3) symmetry applied to the decays D* ZI 7 and D* —> Dtt [17]. Hadronic models tend to 
give values of D* ^ and D* ^ Dtt consistent with (3 near 1 GeV“^ [18].) We will see that the 
contribution in eq. (3.20) is numerically small compared to AMqq computed in the next subsection, 
and the uncertainty in the value of (3 will not be very important for our Hnal results. This issue will 
be analyzed in greater detail in a subsequent paper [ 8 ]. 

We now turn to the contact contribution in eq. (3.4). The ultraviolet sum rule eq. (3.10) allows 
us to parameterize the contact contribution as 




A„ 




(<?! 


mp )2 


(g 2 


E 


IM 


+ 


A 


P'P' 


^p') 


(4 




(3.21) 


where the H’s are linear combinations of the coefficients C and D defined in eq. (3.5). We have 
normalized the H’s so that they are ~ 1 if we identify mp with the hadronic scale that appears in 
these matrix elements. Imposing the infrared sum rule eq. (3.13), we obtain 


App — X App' X H 


pp' 


ip'p'. 


(3.22) 
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It is now straightforward to compute the contact contribution to the electromagnetic mass difference. 
We obtain 


A n^contact 


^ -{Ql - a 


2^ amp ^2 

Stt 


A„ 


1 — + \nx^ 


1 — 


+ Apt pi{\ — x^) 


= -0.02 App> + 0.02 Ap-p- MeV. 


(3.23) 

(3.24) 


We see that for any reasonable value of the A’s, this contribution is numerically negligible, and AMgg 
is given by eq. (3.20) to a good approximation. 

Heavy-quark symmetry equates the electromagnetic mass difference for the P to that of the P*. 
We have (in an obvious notation) 

AM*g = AMgg. (3.25) 


Before continuing the discussion, we note that our final result is rather similar to the prescription 
given long ago for computing hadron electromagnetic mass differences [1], which was essentially to 
evaluate the electromagnetic self-energy with the photon vertices replaced by momentum-dependent 
form factors. However, we emphasize that this prescription was never justified in any satisfactory way. 
In particular, the form factors appearing used in this prescription must be evaluated off shell in order 
to compute the electromagnetic mass difference. Clearly, there are infinitely many ways to continue 
a form factor off shell, and these will give different final results. (In our formalism, this ambiguity is 
contained in the term (!?„) defined in eq. (3.5).) Also, neglecting the contact contributions (contained 
in C in our formalism) was never justified. 

The use of the large-IV limit and the truncation of intermediate states may seem like rather drastic 
approximations. However, we note that precisely analogous approximations are used in the classic 
result [2] for the mass difference: the mass difference is written in terms of vacuum 

correlation functions, which are then saturated with the appropriate lowest-lying 1-particle states in 
the large-limit (in this case, the p and the ai vector mesons). In this way, one can obtain a formula 
similar to the ones obtained here for the heavy meson electromagnetic mass differences [2] [6]: 




3a ml 
1 — mlfm"^ 


9 

Tfl 

In—^ ~ 5.9 MeV. 

ml 


(3.26) 


This is within 30% of the experimental value 4.6 MeV. We therefore adopt this as an estimate of 
the size of the error from the combined approximation of the large-IV limit and the truncation of 
intermediate states. 

As stated in the introduction, the discussion that follows consists mainly of repeated applications 
of the principles described above to compute the remaining contributions to AM. The weary reader 
interested mainly in the bottom line is invited to skip to section 5, which summarizes the main results 
and gives our conclusions. 


3.2. AMqg 

The computation of AMgq can be carried out following the same arguments given in subsection 
3.1, so we will be brief. We can write 


ATQg{v\q^) = -2Y, 


A„W„(g 2 ) 


- A'i+ zO-l- 


- Dn{v ,q ) 


2C{u\q^) 


(3.27) 
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where the functions Wn, Dn, and C are defined as in eq. (3.4), but with Jq replacing J^g. 

The main difference between ATq^ and AT^j^ arises because the form factors of Jq are constants 
in the heavy-quark limit [19]: 


(P(p)|J^(0)|P(p + <z)) = QQU^ (3.28) 

{P{p)\J^mP*{p + q,e)) = 0. (3.29) 


The fact that the form factors are constants even for large momenta can also be understood from 
the point of view of the constituent counting rules [16], since graphs in which the “hard” momentum 
flows directly into the heavy-quark line do not fall off at large \q^\ (see the appendix). From the point 
of view of the large-iV meson theory, we can understand this in the following way: the intermediate 
meson states that contribute to the form factor for Jq in the large-A^ limit are QQ states, and so 
have mass of order 2mQ. These states are integrated out of the effective theory , and their effects are 
correctly taken into account by local interactions. The form factors are therefore pure polynomials in 
q"^ and v in the heavy-quark limit. Because we do not expect the form factors to grow, they must be 
constants, and the constants are fixed by the value of the form factors at g = 0. (Similar arguments 
are made in ref. [20].) 

From eq. (3.28), we then have 


S 

{Dn){q%) = ^ 


a 


(0) 


+ pp{<ll) 


4 + 'ml, 

n(0) 

Te + mis 


(3.30) 

(3.31) 


where the sum over s is over isospin-1 vector mesons that appear in the propagators in figs. 1 and 2. 
As before, and Dn} are constants, while P^^'^ and pP^ are polynomials in g|;. 

We now write ultraviolet and infrared sum rules using the same reasoning used in the previous 
subsection. We have 







(3.32) 

(3.33) 

(3.34) 


As before, we approximate AMg^ by keeping only the contributions of the lowest-lying states; in this 
case, only the P heavy meson state and the p and p' vector mesons are required. 

It is not hard to see that the contact contribution is forced to vanish identically when the sum 
rules are saturated by p and p' vector mesons. The form factor for the light-quark current is given in 
eq. (3.16) in the approximation we are making. Substituting this and the heavy-quark current form 
factor in eq. (3.28) into an expression for AMg^ analogous to eq. (3.5), we obtain 

AMq, cs -Qq{Qu - Qd)^ (3.35) 

^ I -2.5 MeV for Q = c, . . 

1 -bl.2 MeV for Q = b. ^ ' 
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We see that this contribution numerically dominates the light-quark current contribution of eq. (3.20). 
Once again, heavy-quark symmetry gives 


= AMq,. (3.37) 

4. IjiriQ Corrections 

In this section, we consider the 1 / mg corrections to the results obtained above, still working in 
the large-A limit. For mesons containing a c quark, these corrections are expected to be substantial, 
and so it is important to estimate them. 


4.1. AMgg 


To include 0(l/mQ) effects, we must make two types of changes to the formulas of the previous 
section. First, the meson vertices and the masses change by 0(l/mQ) effects. Second, the form of 
the heavy-meson propagators is modihed. This modification can be thought of as a recoil correction: 
it is completely kinematical in origin, and hence determined with no dynamical input. Formally, 
this can be thought of as writing a non-relativistic expansion of the fully relativistic propagator or 
(more correctly) an expression of reparameterization invariance [11]. The result is that eq. (3.5) in 
the previous section is replaced with 


AMgg = 


f 1 

J (27r)4 g| 



(A„ -b (5An(gg))VF„(-g|) 
(<?£4)^ + Al 


{Dn){ql 


+ {c){q%) 


(4.1) 


where 


SAniqE) 


q% + Al f 2Al \ 
2M { qh + AlJ- 


(4.2) 


The functions (C) and {Dn) are given by formulas identical to eqs. (3.7) and (3.8), except that the 
coefficients are now given in a power series in I/itiq: 


Cif = + Oil/ml), 

AMgg = AM® + AM(i) + 0(l/m|), 

etc., where AMgg\ ... parameterize the 0(l/mQ) corrections. 

We obtain ultraviolet sum rules by imposing the condition that the integral for AMgg converges 
order by order in IjiriQ. Because the form factors fall off as ~ 1/q^ (see the appendix), the form-factor 
contribution for each n converges by itself, and the sum rules again constrain the contact contribution 
alone: 

E (dl) + E E = 0, (4.4) 

r,s n r,s 

E + E E d = 0,1. (4.5) 

r,s n r,s 


To get the infrared sum rules, we note that the recoil corrections do not affect the leading behavior 
of the integrand as qs —*■ 0. However, there is a new feature that enters at 0(1 /toq) because 


AT{p,q)^{Ql-Ql) 


2TriS{qo) 


4 

M 


as q^O. 


(4.6) 
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This results by writing a low-energy effective theory where only the P and the P* appear. In this 
theory, the “extra” 1/M term arises from the term 


= -Lpt£) 2 p 

2M 


(4.7) 


where the coefficient is fixed exactly by reparameterization invariance. (Other interactions that arise 
at order l/mq vanish as g—>0.) Because of this, the infrared sum rules are 


E 


Q 


(0.0) 


E 


a 


r,s 

(0,1) 




EE 


D 


(0.0) 




= 0 . 


EE 


D 


(0.1) 


= -{Ql-Ql) 


(4.8) 


(4.9) 


00r'"-ls „ r,s '""Or'"-Is 

We have once again used the fact that the form factor contribution coming from Wniq"^) already has 
the correct infrared behavior, except for the term in eq. (4.7). 

These sum rules can be saturated by the same set of states as in the previous section: p, ui, p', 
Lo', P, and P*. The form factor contribution is convergent by itself, and is given by 

12 n2\^,^2 

form factor ^ 

‘-qq 


^^(1) form factor _ _ (Qa " 1 + 2x^ In 

47rM (1 — x^)^ 


3(Q^ - Ql)al3^ml 


Stt 


{M* - M) 


1 -I- In 

(1-:e2)3 


(4.10) 


nip 2{1 — x^)-\-{1-\-x^)\nx 


21 


2M (l-a;2)3 

Using the infrared sum rule, the contact contribution can be written (compare to eq. (3.21)) 

^ Ay (1) contact) (^2^) ^ iQu~ Qd)''^p 


M {q^ + 772,2)2 

{Ql - Ql)ml 


A 


( 1 ) 

pp 


^PP' 


.( 1 ) 

^P'P' 


(g|-fTO 2)2 (q| +to 2 )(^| _^^2 ) (^|^_^ 2)2 


(4.11) 

The constants ~ irip/M parameterize the 0{l/mQ) corrections to the unknown contact contri¬ 
butions defined in eq. (3.21). The ultraviolet sum rule then gives 


3(1)__ MaA) 


PP 


P P 


and we obtain 


A A^(l) contact 


{Ql - QDarrip 
8ttM 
[Ql - Ql)am^ 


Stt 


P 2 
X 


' (i)l-a :2 + lna :2 ' 

^pp' +^p’p’A X) 


(4.12) 


(4.13) 


We find that the corrections that depend on the r4’s are numerically negligible for reasonable values 
of AA') . Our result is therefore 


am^a 


0.088-f 0.028 
0.031 -f 0.0094 


/? 


1 GeV^ 
(3 

1 GeV^ 


MeV for Q = c, 


MeV for Q = b. 


(4.14) 


13 



Comparing to the lowest-order results in eq. (3.20), we see that the l/mg corrections are 35% for the 
c system and 10% for the b system for /3 = 1 GeV“^. 


At order IjniQ the P* electromagnetic mass difference is no longer equal to that of the P. 
Applying the same method, we obtain the isospin-violating hyperfine splitting due to the light quark 
charges. 


am;, - AM,, ^ 


{Ql - Ql)a(3^m- 


27T 

[Ql - Ql)af3m- 


6 


(M* - M) 

hp'-/3) 


1 -I- 2x'^ In — x'^ 
(1-X2)3 


(1 -f xY 


(4.15) 


where f3' measures the strength of a P*-P*-^ coupling. Heavy-quark symmetry gives (3' = (3 + 
0(l/mg), but we do not have any experimental information about the difference (3' — (3- We therefore 
obtain 


-0.16-0.015 


am;, - AM,, 


/3 


1 GeV" 


0.3 P 


-0.054-0.0052 


P 


1 GeV 


P'-P 

0.1 p 


MeV for <5 = c, 
MeV for Q = b, 


(4.16) 


where we have normalized the value of /?' — /? to a representative value. Comparing to eq. (3.20), we 
see that the l/mg correction to AM;, is 20% for the c system and 5% for the b system for the values 
of P and (3' used to normalize these expressions. 


4.2. AMQg 


When we include the 0(l/TOg) corrections, AMqq is given by a formula just like eq. (4.1), 
but with Jq replacing J^g. The functions (C) and (Dn) are given by formulas like eqs. (3.30) and 
(3.31), but with the coefficients expanded in powers of l/mg, as in the last section (see eq. (4.3)). 
Using reasoning similar to that of the previous subsection, we obtain ultraviolet sum rules of the 
same form as eqs. (4.4) and (4.5), and infrared sum rules of the same form as eqs. (4.8) and (4.9). 
The contact contribution is uniquely determined when the sum rules are imposed. Performing the 
necessary computations, we obtain 

QQiQu-Qd)aml \nx^ ^ Qq{Qu - QYaPPqml 1 
~ 2ttM 1-x^ 2 :e(1 + x)’ ^ ’ 


where Pq is the coupling of the heavy-quark current to P and P*, normalized like P in eq. (3.17). In 
the heavy-quark limit, 

pQ=^[l + 0{aYM))]. (4.18) 

The corrections are expected to be of order as(mg)/7r, which is 0.15 for the c quark and 0.08 for the 
b quark. Neglecting these corrections, we obtain 




-0.44 -k 0.73 (- - —i- ) MeV for Q = c, 

VI GeV^ J 

+0.077- 0.13 1 —^ ) MeV for Q = 6. 

Vl GeV^ J 


(4.19) 


Comparing to eq. (3.36), we see that this is a 10% correction for the c system and a 15% correction 
for the b system for /3 = 1 GeV“^. (The small size of the corrections for the c system is a result of 
cancelations that depend on the value of p.) 
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Similarly, we can compute the P* isospin splittings: 


- AMq, ~ 


2Qq(Q„ - Qd)aP(3Qm^p 1 


which gives (again using eq. (4.18)) 


AM^, - AMq, 


-0.99 


+0.17 


P 


1 GeV 


P 


1 GeV 


x{l + x) ’ 


MeV for Q = c, 


MeV for Q = b. 


(4.20) 


(4.21) 


5. Summary and Conclusions 


We now summarize our results. We have computed the heavy-meson electromagnetic mass dif¬ 
ferences by working in the large-A limit, where the electromagnetic mass differences are given by 
a convergent sum over an infinite number of 1-particle intermediate states. We obtained sum rules 
relating the matrix elements appearing in this sum that enforce the correct ultraviolet and infrared 
behavior on the electromagnetic amplitudes that appear in the calculation. All of this is a rigorous 
consequence of QCD (in the large-A limit); however, in order to get numerical results, we truncated 
the infinite sum by keeping the smallest number of intermediate states that are capable of giving a 
consistent description of the matrix elements which appear in the sum. We argued that these approx¬ 
imations are similar to the ones made in the classic calculation of the 7r“'"-7r*^ mass difference, which 
works to 30%. Making these approximations, we find that the numerically dominant contribution 
comes from the heavy-quark current, and we obtain (see eqs. (3.20), (3.36), (4.14), and (4.19)) 


{Mdo - ^ 


(Mb+ - Mbo)^“ ~ 


-2.4 + 0.74 


+ 0(l/m^). 


+ 1.7-0.13 


+ 0(l/mg). 


P 


1 GeV^ 


P 


1 GeV^ 


- 0.012 


- 0.03 


P 


1 GeV^ 


2n 


MeV 


P 


1 GeV 


-1 


MeV 


(5.1) 


(5.2) 


Here, /3 is a matrix element that measures the strength of the P*-P-^ coupling (see eqs. (3.15) and 
(3.17)). The coefficients of the terms linear in /? are 0(1 /toq) and have perturbative QCD corrections 
of order q;s(toq) (see eq. (4.18)). Similar results are also obtained for the vector mesons. Heavy quark 
symmetry gives AM* = AM + 0(1 /toq), and we obtain (see eqs. (4.16) and (4.21)) 


[MjjtO — Md->+) ~ {Mjyo — Md+) 


EM 


-0.16-0.99 

+ Oil/ml) 

{Mb^+ - Mb-o)™ ~ {Mb+ - Mbo)™ 


P 


+ 


-0.054 + 0.17 

Oil/ml), 


1 GeV 


P 


-1 


- 0.015 


P 


1 GeV 


-1 


P'-P 

0.3 P 


MeV 


(5.3) 


1 GeV 


-1 


- 0.0052 


P 


1 GeV 


-1 


P'-P 

0.1 p 


MeV 


(5.4) 


where P' = P + 0{l/mQ) is another unmeasured matrix element (see the discussion above eq. (4.16)). 
As above, the coefficients of the terms linear in /3 are 0(l/mQ) and have corrections of order asirng). 
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For (3 near 1 GeV~^, the heavy-quark expansion appears to be working very well for the b system, 
and moderately well for the c system. We believe that this is sufficiently encouraging to consider the 
phenomenology of these results in detail in a subsequent paper [8]. For the present, we hope that 
the systematic approach taken in this paper is the starting point for further progress on this classic 
problem in hadronic physics. 
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Appendix A. Ultraviolet Behavior of Heavy-Meson Form Factors 

In this appendix, we consider the ultraviolet behavior of matrix elements of the form 


{P{p)\Jq{0)Hp + q)), (A.i) 

where n is a heavy-light meson. We work in the the l/mq expansion. This means that we consider 
spacelike momenta q in the limit 

Aqcd ^ (A.2) 

The ultraviolet behavior of these form factors can be determined from the constituent counting rules 
[16]. The idea is that for large spacelike momentum transfer, the leading behavior of the form factor can 
be obtained by power counting the hard momentum flow through constituent Feynman diagrams with 
appropriate kinematics for the initial and final state constituents. For the case under consideration, 
the leading graphs come from a single hard gluon exchange, as shown in the first two graphs in fig. 3. 
Evaluated in the heavy-quark efi'ective theory, these graphs give a contribution of the form 


{k + qY{xqY 


•(/>'( fc') 7^ (|(! — x4)i/i'tl}{k) 

{k — xqY{xqY 


where x is the fraction of the hard momentum q flowing into the heavy quark. The average is taken 
over k and x. The dominant contribution will come from regions of integration where almost all of the 
hard momentum flows into the heavy quark (i.e., k ~ Aqcd and a: = 1 -I- 0(Aqcd/wq))- The reason 
is simply that if hard momentum flows into the light quark, the final state of the hard scattering 
subprocess consists of the heavy quark at rest and a light quark with momentum much larger than 
Aqcd; such a state is expected to have very small overlap with a heavy meson at rest. Therefore, 

fig. 3 yYl'" - '<Pik)) +-^{W{k') ■fpik)} H-, (A.4) 

and we have 


{P{p)\J^mP{p + q))r 

Aqcd 

9^ ’ 

(A.5) 

{P{p)\J^mP*{p + q))^ 

Aqcd 

9^ ’ 

(A.6) 


etc. For the P elastic form factor, we have used the fact that rotational symmetry implies that the 
first term in eq. (A.4) is proportional to q^v^v'' = 0(l/mQ), since the fact that the initial and final 
states have the same mass forces q ■ v = —q^/2M. In general, it is clear that matrix elements of the 
form of eq. (A.I) fall off at least as fast as ~ l/?^. This is the result quoted in the main text. 

A complete analysis of the form factors would include a discussion of higher-order graphs such 
as the last one in fig. 3. We will not give a detailed analysis of this issue here. However, we expect 
graphs such as this to factorize, so that regions of integration where the gluon momenta are soft can 
be absorbed into “wavefunction” corrections, while the remaining hard contributions have the same 
asymptotic behavior as the contributions analyzed above (up to logarithms). For an example of this 
type of analysis, see ref. [21]. 
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Figure Captions 


Fig. 1. Contributions to the Compton amplitude T in the large-iV limit. The sum over n runs over 
excited heavy mesons, while the sum over r and s runs over vector mesons. The shaded blobs on the 
graphs involving a sum over n are heavy meson form factors; see fig. 2. 


Fig. 2. Contributions to heavy-meson form factors in the large-A'^ limit. The sum over r runs over 
vector meson states. 


Fig. 3. Contributions to the quark scattering amplitude used to determine the asymptotic behavior of 
the light-quark current form factors. The thick line is the heavy quark, and the curly line represents 
the gluon propagator. 
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This figure "fig 1-1.png" is available in "png" format from: 


http://arXiv.org/ps/hep-ph/9502259vl 
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